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The electronic Lorentz theory is employed to explain the optical properties of planar split-ring metamaterials. 
Starting from the dynamics of individual free carriers, the electromagnetic response of an individual split-ring 
meta-atom is determined, and the effective permittivity tensor of the metamaterial is calculated for normal in- 
cidence of light. Whenever the split ring lacks in-plane mirror symmetry, the corresponding permittivity tensor 
has a crystallographic structure of an elliptically dichroic medium, and the metamaterial exhibits optical prop- 
erties of planar chiral structures. Its transmission spectra are different for right-handed vs. left-handed circular 
polarization of the incident wave, so the structure changes its transmittance when the direction of incidence is 
reversed. The magnitude of this change is shown to be related to the geometric parameters of the split ring. The 
proposed approach can be generalized to a wide variety of metal-dielectric metamaterial geometries. 
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I. INTRODUCTION 

The concept of metamaterials have drawn a strong inter- 
est ever since fabrication and characterization of such struc- 
tures became feasible. The ability to engineer the metama- 
terial elements ("meta-atoms") in a largely arbitrary fashion 
adds a whole new dimension of freedom in material science. 
Artificial materials show promise for a wide range of un- 
usual physical phenomena, some of which are rare or absent 
in nature. Examples include negative refraction,^ hyperbolic 
dispersion relation resulting in anomalously high density of 
states in a wide frequency range,- or support for transforma- 
tion opticsr-"- offering superior degree of control over light 
propagation and guiding. 

More specifically, if the shape of the meta-atoms is chi- 
ral, i.e., when the meta-atom cannot be superimposed with 
its own mirror image, such metamaterials resemble naturally 
occurring optically active media and outperform them by or- 
ders of magnitude. More recently, following a seminal pa- 
per by Plum et al,^^ planar chiral metamaterials (PCMs) 
were introduced. They consist of planar meta-atoms on a 
flat substratei^— that cannot be superimposed with their in- 
plane mirror image without being lifted off the plane. In other 
words, truly chiral metamaterials possess distinct 3D enan- 
tiomers while PCMs do not, but can be said to possess distinct 
2D enantiomers [see Fig. [Tib)]. 

A notable feature of the PCMs is their asymmetry in prop- 
agation of electromagnetic waves incident from opposite di- 
rections or having different (right- or left-handed) circular 
polarization.^^ This directional asymmetry results from the 
fact that the wave propagating in the opposite direction ef- 
fectively interacts with the mirror image of the original struc- 
ture. Therefore it is a manifestation of asymmetry between 
2D enantiomers. This makes PCMs distinct from 3D chiral 
metamaterials^'^' where 3D enantiomers do not flip when 



the wave propagation direction is reversed. 

Directional asymmetry is also known to exist in other en- 
vironments such as the Faraday cell. In this case, however, it 
is usually attributed to non-reciprocity due to the presence of 
magnetic field, while PCMs do not have such non-reciprocity. 
Further analysis reveals that both in Faraday media and in op- 
tically active liquids the waves whose polarization state does 
not change during propagation (polarization eigenstates) are 
circularly polarized and counter-rotating, i.e., they come in 
pairs of eigenwaves one of which has right-handed (RH) and 
the other left-handed (LH) circular polarization. In contrast, 
PCMs have polarization eigenstates that are elliptical and co- 
rotating}^ meaning that their handedness is the same for both 
eigenwaves in the pair. 

Thus, while 3D chiral metamaterials are simply an analogy 
of bi-isotropic (e.g., optically active liquids) or bianisotropic 
(gyrotropic) media, PCMs clearly represent a distinct type of 
electromagnetic materials, apparently without an immediate 
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FIG. 1: (Color online) (a) Schematic of a bilayer (3D) vs. planar (2D) 
chiral metamaterial; (b) illustration of 3D and 2D enantiomers 
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naturally occurring counterpart. Hence, explaining the exotic 
optical properties of PCMs certainly poses an exciting prob- 
lem in theoretical electrodynamics. On an abstract crystallo- 
graphic level, it has been shown recently that a combination of 
birefringence and circular dichroism provides a basis for op- 
tical manifestation of planar chirality, and that an elliptically 
dichroic medium should exhibit characteristic optical proper- 
ties of a PCM. 14 

On a more involved, microscopic level, several explana- 
tions of PCM operation have been developed. Most works 
attribute the dichroism to polarization- sensitive excitation of 
dark plasmonic resonant modes in meta- atoms. In a vari- 
ety of PCM designs, this polarization sensitivity is attributed 
to magnetic dipole excitations that become "trapped" due to 
poor coupling with the incident waveJ^ai^^i^ Other accountsii 
build up on earlier multipole treatment of metamaterialsi^ to 
express the meta-atom's dichroic response based on collective 
dipole oscillations in its individual segments. Still other works 
attempt to formulate the theory of PCMs by treating the meta- 
atom as a plasmonic oligomer and calculating its response in 
coupled dipole approach. ^^-^^ 

However, an ab initio theoretical description of PCMs, 
which would explain their dichroic behavior on a microscopic 
level and independently of a particular PCM geometry, is still 
missing. A first attempt, based on considering the response of 
the meta-atom in terms of the dynamics of individual electrons 
(Lorentz theory), was made in Ref.[3. That attempt, however, 
was only successful in arriving at the right form of the meta- 
material's effective dielectric permittivity tensor, reproducing 
qualitatively correct optical properties of a PCM. A more rig- 
orous, quantitative description with a detailed analysis of the 
applicability range for the Lorentz theory is still called for. 

In this paper, we make a further step towards an ab ini- 
tio description of metamaterials and present a detailed route 
to arrive at the chiral properties for the asymmetric split-ring 
PCMs.—^ We start by considering the dynamics of free elec- 
trons in an arc- shaped metallic element in presence of an ex- 
ternal electromagnetic wave. The electron motion is shown to 
be determined by (i) screening forces caused by the charge re- 
distribution in the external field, and (ii) electromotive forces 
caused by currents flowing in the segment as the charges re- 
distribute. These forces exert spring-like and inertia-like ef- 
fects on an electron, respectively, and can be associated with 
capacitance-like and inductance-like contribution in an equiv- 
alent LC-circuit. Together, these two kinds of forces give rise 
to a particle plasmon resonance. Its properties are derived di- 
rectly from geometric dimensions of the arc segment and ma- 
terial properties of the metal without the need for phenomeno- 
logical parameters. By considering two such segments in 
close proximity, the electromagnetic response of a split-ring 
meta-atom is evaluated. 

In order to connect the meta-atom's response to the exotic 
properties of PCMs described previously, a standard homog- 
enization procedure^i^^^ is then employed to derive the ef- 
fective dielectric permittivity tensor for the PCM. Since the 
PCM in question is a surface rather than a bulk materialr^i^i 
this procedure should not be regarded as a true homoge- 
nization, and a multipole-expansion approach should give a 



more accurate physical description i ^1^5^26 ^^{\\^ interpreting 
a metafilm as an effective medium can be feasible^ and the 
resulting effective permittivity tensor is shown to be valid for 
normal incidence of light in the frequency range of the fun- 
damental resonance of the split ring. Moreover, it turns out 
to have the crystallographic structure of elliptically dichroic 
media, as would be expected for planar chiral materials. A 
similar approach was recently applied to calculate the THz 
field enhancement in a nanoslit,^^ the results showing a good 
agreement with experiment. 

The optical properties (transmission and reflection spectra) 
of the PCM are then calculated using wave operator based 
extension of standard transfer-matrix method^^^^i In agree- 
ment with previous theoretical and experimental results r^^i^ 
the spectra are sensitive to whether the incident wave has left- 
or right-handed circular polarization. This difference, which 
translates to directional asymmetry, can be used to quantify 
the strength of planar chiral properties. The dependence of 
this strength on the geometrical parameters of the split ring is 
analyzed. When the split ring is symmetric, i.e., when there 
are no distinguishable 2D enantiomers, optical manifestation 
of planar chirality is seen to vanish as required by symmetry 
and reciprocity constraints.^^ Maximum chiral properties are 
observed when the enantiomers are most distinct. 

The proposed approach can be relatively straightforwardly 
extended to a wider variety of shapes for planar meta- atoms. 
Moreover, the model lends itself to an extension along the 
lines ofRefs.llTHTsI to the cases of oblique incidence and non- 
planar meta-atom shapes. The results obtained can also be 
generalized from a single PCM to PCM-based multilayers, 
inasmuch as such generalization can be performed.^^ 

The rest of the paper is structured as follows. In Section |II1 
the Lorentz theory is employed to arrive at the electromagnetic 
response of a single split-ring meta-atom. Section Hill follows 
with the procedure to arrive at the effective permittivity tensor 
of the PCM. The structural properties of this tensor are ana- 
lyzed, and calculation of the PCM's optical spectra is outlined. 
In Section |IVl the results on calculated spectral properties of 
split-ring PCMs are presented and compared with numerical 
simulations. In Section |Vl the relations between chiral prop- 
erties and meta-atom geometry are systematically analyzed. 
Finally, Section |Vl] summarizes the paper and outlines future 
extensions for the proposed theory. 



II. RESPONSE OF A SPLIT-RING META-ATOM 

As an example of a planar chiral metamaterial, we consider 
a two-dimensional array of chiral split rings (CSRs). The 
corresponding meta-atom is shown in Fig. [S] It comprises 
a metallic ring broken into two segments in an asymmetric 
fashion, so that two 2D enantiomers can be distinguished [see 
Fig- Eh)]. CSR metamaterial was chosen for its relative geo- 
metrical simplicity and for availability of previous experimen- 
tal results. 1^ 

The lateral width of the ring d is assumed to be much 
smaller than its radius R, and its thickness h is even smaller. 
The ring sits atop a thicker dielectric substrate. The metal of 
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the ring is taken to be copper with complex dielectric permit- 



tivity £m,^ and the permittivity of the substrate is e^. Both 
materials are non-magnetic. A monochromatic incident wave 
with electric field E = Eq exp(— ico^) is assumed to illuminate 
the material, (0 being the angular frequency of the wave. 



A. Dynamics of the electrons 

We begin by considering the motion of an electron in a 
finite- sized metallic inclusion. Each electron is affected by the 
driving force eE originating from the external electric field, 
the friction force — myv due to losses in the metal, the screen- 
ing force Fc resulting from other electrons as they are redis- 
tributed in the inclusion by the external field, and the electro- 
motive force ¥l due to currents produced by electrons as they 
move under the action of the field. Here, m and e <0 are elec- 
tron mass and charge, v is its velocity, and 7 is the decay fre- 
quency for the metal (a friction coefficient for electrons). By 
determining the dynamics of the electrons in the metal and by 
averaging this dynamics over the meta-atom, its polarizability 
can be determined, and one can introduce effective dielectric 
parameters of the entire metamaterial. 

It should be noted that averaging and homogenization are 
two distinct procedures. From the averaging of the micro- 
scopic parameters we get the average displacement of elec- 
trons and the overall polarizability of the metal ring segments. 
Specifically, the averaging yields the resonance frequency and 
spectral shape. Once the polarizability of the metal segments 
is known, the effective medium parameters can be introduced 
by homogenization. To be able to connect the properties of 
PCMs with those of bulk planar chiral media, we employ the 
standard approaches used for bulk metamaterials. However, it 
should be realized that the resulting effective parameters can 
be attributed to the PCMs under study only in specific cases 
(for the normal incidence of light). 

The screening force Fc is expected to depend heavily on 
the geometry of the metallic inclusion. Under the action of 
the field, all the electrons are displaced, and uncompensated 
charges are accumulated at the edges of the metal in the di- 
rection of the field. If the distance between the edges in that 
direction is small and the edges are wide enough, the induced 
charges produce the field similar to that inside a capacitor, i.e. 
it is mostly homogeneous and compensates the external field. 
Therefore it can be assumed that the electrons simply won't 
move between such edges. So, for our geometry we can ne- 
glect out-of-plane and radial electron motion as long as d and 
h are small. 

The displacement of an individual electron therefore de- 
pends on time and on the position in the ring: ^ = ^(^, 9). It is 
subject to the equation of azimuthal motion for the electrons 

(C = 5CM): 



m'C + myt, = • (eE + Fc + Fz,) . 
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FIG. 2: (Color online) (a) Top view and (b) side view of a CSR unit 
cell; (c) schematics of the induced charge redistribution in a CSR: 
A^i 2 = ^1,2^^' where A»S = hd. The parameters are as in Ref.[T2l 
(cell size: L = 15mm, ring radius: R = 6 mm, ring width: J = 0.8 
mm, ring thickness: h = 35 /im, substrate thickness: D = 1.6 mm, 
substrate material: dielectric with =4). The angles determining 
the CSR composition are ai = 140°, 0^2 = 160°, Pi = 40°, = 20°. 



and charges at the tips of the ring segments. The charges cre- 
ate a "restoring force" Fc like the one in a harmonic oscillator 
(e.g., by Hooke's law). The electromotive force Fl caused by 
the current is the reason of the additional non-dynamic iner- 
tia of the electrons. As in the ordinary classical mechanics, a 
particle under these forces will be subject to resonant motion, 
which should manifest itself as resonant properties of the ef- 
fective medium. 

Usually, the electron displacement is small and in good ac- 
cordance with the external field, so for a sub wavelength- sized 
inclusion it can be considered as position-independent. In this 
case, = ^{^)- However, this behavior is expected to 

break down in the vicinity of a particle plasmon resonance 
where the electron shift can be substantial and more in accor- 
dance with the resonant mode. Therefore, ^ becomes very 
sensitive to the (p-dependence. At the fundamental resonance, 
it is reasonable to assume the standing-wave positional depen- 
dence in the two ring segments 



Ci(^ (p) = si (t){n/2) sin(;r((p - (pi)/ai), 
C2(^ (p) = S2(t){n/2) sin(7r((p - 93)70^2), 



(2) 



Equation of motion introduces resonant features in the re- 
sponse of the structure. Electrons are pushed by the external 
driving forces, which induce both electric currents in the ring 



where ^1 = (l/ai)/^^^ ^^j^ ^^d ^2 = (1/0^2) J^^' ^2^9 have 
the meaning of an averaged electron displacement in the cor- 
responding segment. It is seen that the electrons do not move 
near the tips of the split rings as 9?) =0' ^ = 1 5 • • • ? 4 (see 
Fig. [21), and maximally shift in the middle of the segments. 
While the positional dependence in can be safely neglected 
in the off-resonant situations, it will be shown below that it is 
required to correctly predict the resonant frequency of a CSR. 

The effective response of the entire meta-atom is obtained 
by averaging over all available azimuthal angles where metal 
is present, i.e., for (p from (p\ to (p2 and from 93 to 94 [see 
Fig. [2tc)]. With the brackets (A) denoting the averaging of 
A over 9, the equation of motion for the averaged electron 
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displacement s takes the form 

ms^mrs = e{{e^ • E^)) + {e^ • Fc) + {e^ • Fl) . (3) 

We assume that the metal is fully embedded into the substrate 
medium, so external field E coincides with the field in the 
ambient medium E'^. 



B. Screening force and induced charges 

Since the accumulation of charges along the ring edges only 
serves to prevent the electron motion in the radial direction, 
the screening force Fc inside the ring segment can be ap- 
proximately described as the Coulomb force generated by the 
charges at the tips of the segment, labeled by the points (pi and 
(p2 ((P3 and 94). This assumption (which is valid outside the 
immediate vicinity of the tips) is further substantiated by the 
fact that charges do accumulate mostly near the tips of metal- 
lic objects. In a CSR, we need to take into account the tips of 
the other segment in the same split ring [see Fig. [SJc)], as well 
as the influence of the neighboring CSRs. 

The screening force on electrons in a ring segment from its 
own tips is essentially similar to that in a metallic rod with 
width d and thickness h. For the rod placed in dielectric with 
e = ej, the force is equal to 

F^^P(x) = -l^ arctan f^J£/^=y (4) 

where Z is the surface charge density at the tip facet of the rod 
with area AS = dxh , and x is the distance from the center of 
the facet to an observation point. 

To avoid unphysical behavior of the force F^^^ very close to 
the tip, the positional dependence of the displacement ^(9) in 
Eq. (O has to be taken into account as L = eN^{t^ (p). 

The expression given by Eq. ^ is approximately valid if x 
is small enough so that the arc segment with length x is not 
significantly different in shape from a rod with dimensions 
X X d X h. It can be assumed so if the line connecting the 
observation point with the tip lies wholly within the metal, 
i.e., if < X < 2VRd for /z <C J. 

Outside of that range, Eq. (|4]) no longer holds. However, for 
the observation point far away from the tips one can regard the 
charges accumulated at the tip as one point charge and calcu- 
late the screening force according to the Coulomb law. Unlike 



F^^P, the forces "in the center" of an arc segment (labeled F^^^) 
are not dominated by contribution from any particular tip, so 
both tips of the segment in question, as well as the tips of 
the neighboring segment in the ring and from the neighboring 
rings need to be taken into account. As a result, we get the 
following expression for the upper segment 

Fr = Fc-e^ 

\Rer{(p)-hnm-Rer{(p2)\^ J 
_^ eZ2ASy f ecp{(p){hnm^R^ri(p3)) 

^d \R^r{(p) - Km - Rer{(p3)\^ 

e(p{(p){Km^R^r{(p4)) \ 
\Rer{(p)-Km-R^r{(P4)\y' 

where Km = n^x + is the position of a CSR in the (x, 
y) plane, characterized by a couple of lattice indices n and 
m. The charge densities for the upper and lower segments are 
Zi = eNC^i (t, (p) and Z2 = eNt^2{t, 9), respectively. 

To obtain the averaged screening force for the meta-atom, 
we assume that it is given by E q. (|4| for < x < 2\^Rd (i.e., 
within an angle 5(p = l^djR from each tip, see Fig.[2lc)), 
and by Eq. (|5]) elsewhere. Accordingly, the averaging over the 
upper segment results in the formula 

(e^Fc)i = — / F'^^dcp + / Frdcp^ 

r F^pjcpY (6) 

J(p2-d(p J 

Similar expressions can be obtained for the lower segment. 
These integrals cannot be evaluated in closed form, but can be 
easily calculated numerically. The resulting screening force 
for the two segments in the ring can be finally expressed as 

^1,2 = — ^ll,22'5'l,2 -^12,21'5'2,l, (7) 

where kij are coefficients obtained from the integrals in 
Eq. ©. They have the meaning of "stiffness" coefficients in a 
mechanical oscillator, and are given by 
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C. Electromotive force 



where 



The force = eEem appears because the electron motion 
along the split ring under the action of the external electric 
field can be regarded as currents in the loop and induces sec- 
ondary magnetic field, which in turn penetrates the closed 
contour of the loop and creates an electromotive force. This 
force is directed in opposition to the induced current in accor- 
dance with Lenz's rule. 

To determine F^, we start from the Maxwell equation in 
integral form 



1 



(9) 



The integration is performed along the loop of the ring C, 
while S = ttR^ is its area and is the z-component of mag- 
netic field. Assuming homogeneous electromotive electric 
field over the contour and taking the length of the loop as 271 R 
and its area as 5" = ttR^, the force affecting an electron in a 
CSR is 



2c 



(10) 



Since the magnetic field of a normally incident wave has no 
z-component, can only originate from the current flowing 
along the loop. The field can be estimated from the Biot- 
Savart law 



B(ro) = / 



Id\ X r 



(11) 



where r is the radius-vector from the current element Id\ to 
the observation point fq. The current / in the metallic parts of 
the loop is I\^2 = eN^i^2^S = eN^i^2hd. 

We then assign the magnetic field in the loop to be the mag- 
netic field at the center of the split ring. The magnetic field 
created by the current in each segment is summed: 



1 / 
cR 

so the electromotive force finally equals 

Fl = {ecp ■ Fl) = -^J^ [^I'^i + ^2S2] . 



(12) 



(13) 



The force Fl is proportional to the electron's acceleration 
which allows us to regard the coefficient in front of s\^2 as an 
effective mass. 



D. Equation of motion 

By substituting the screening Eq. d?]) and electromotive Eq. 
([T3]) forces into equation of motion for the electrons © in 
each segment, we derive the following coupled differential 
equations: 

mns'i -^mnh^yh +^ii^i +^i2'^2 = ^((e(p)iE^), 
m2lS\^Jn22S2^7S2^k22S2^k2lSl = e{{e(p)2^^) , (14) 



mil = Jn-\-m 
mi = 



aiCOphd 



ai (Ophd 



m\2 = m 



m22 = fn-\-m 



a2C0phd 
a2(0phd 



(15) 



Here (e^)i = ^^^J^e^dcp and (e^)2 = ^^J^^^^dcp 
denote the averaging of the vector e,^ over the upper and lower 
segment, respectively (see Fig.©. 

For the geometrical parameters of CSRs used, it can be es- 
timated that maiCOphd/^nc^ m, which leads to mn ^ m2i 
and mi2 ^ m22- In other words, the effective mass of an 
electron is dominated by the contribution of the electromo- 
tive force rather than by the dynamic counterpart. However, 
electromotive and dynamic masses can become comparable at 
the nanoscale. 

Note that the effective mass ntij can be regarded as induc- 
tance, while the coefficients kij have the meaning of inverse 
capacitance. In this picture, it is seen that the proposed model 
coincides with a well-known effective-circuit (LC) model 
with the effective inductance and capacitance calculated from 
the material and geometrical properties of the CSR. This re- 
sult is not surprising because the LC model is expected to be 
valid for millimeter- sized resonators at the gigahertz operat- 
ing frequencies. It also follows that the LC-model remains 
applicable as long as AS = hd:^ Snc^/{aiCOp) ~ 10~^^ m^. 

In a CSR, the coupling between the two segments is real- 
ized by means of the effective masses mi2 and m2\ and the 
"stiffness" coefficients kn and k2i. As we have just seen, all 
mass coefficients rriij are almost identical. It is not the case for 
the coefficients ku, which contain the contribution in the vicin- 
ity of the tips F^^P [see Eq. Q] . This contribution is greater 
than the forces in the central part of the ring F^^^ by about 10 
times, so it can be considered dominant. On the contrary, the 
coefficients kij^i contain only F^^^, so kij <C ku. Moreover, 
since ku muCO^ ^ rriijCO^ > kij, so the condition F^^^ <C F^^^ 
allows the contribution of kij^i to be entirely neglected. Thus, 
we can simplify the expressions for the remaining stiffness 
coefficients ku: 



^11,22 




)^ rKd(p 

^ / sin(7rx/(ai,2^)) 

,2^ ^0 

hd/4 



dx. 



(16) 



Note that both tips for each segment are taken into account. 

Eqs. ([T4l) are essentially equations of motion for two cou- 
pled oscillators with a time-harmonic driving force. There- 
fore, we look for a solution of these equations (O in the form 
si 2 = exp(— i()L)r)/i^2, where (0 is the frequency of the incident 
wave and l\ and I2 are constants. Then the equations are sim- 
plified as 

(-miico^ - i(oy^kn)h + {-mnco^ ^ki2)l2 = e{e(p)iE'^, 

{-m2iCO^^k2i)h^{-m22CO^ -icoY^k22)h = e{ecp)2E'^. 

(17) 
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This system is easily solved with respect to /i and I2. So, 
we find the averaged displacements of electrons in both CSR 
segments of the ring: 



/i 



1 



47lNe 



h = ^tei(co)(e^)i +;^22(co)(e^)2)E^ (18) 



where 



Zl2(ffl) ^ 



-mco'^{m22(0^ + iyco - ^22) 
d'o ' 
mC0p{mi2C0^ — ki2) 
d'o ' 
mC0p{m2iC0^ — k2i) 
d'o ' 
-mcOp (mi 1 (0^ + iyco -ku) 
d'o ' 



Do = {mn(0^ + iY(0-kn){m22(0^ + iY(0-k22) 

-(ot12«^ - ki2){m2i(0^ -k2i). 



(19) 



E. Polarizability of the meta-atom 



Finally, to arrive at the effective polarization of the unit 
metamaterial cell, we present it as a sum of the polarization 
of both CSR segments and of the surrounding dielectric: 

47r{F) = {l-pi-p2){£d-l)E' 

+ 47rA^^(/?i(e(p)i^i+/?2(e(p)2^2) , (20) 

where /?i,2 = Rda\ 2h/L^D are filling factors of the two metal- 
lic segments. By substituting the average displacement (fTSl) 
into Eq. ([2Q|) . we get 

47r(P) = {\-p,-p2){ed-m' 

+ /^2(e(p)2 0U2i(e(p)i+;^22(e(p)2)]E^, (21) 

where (E'") is the thickness-averaged electric field in metal, 
E^ is the field in dielectric, and a b denotes an outer (tensor 
or dyadic) product between two vectors. 

The average electric field in the metamaterial cell approxi- 
mately equals the electric field in dielectric, (E) ^ E^. So, the 
polarization of the cell can be expressed as 

4;r(P)=47r;^(E) = [(8^-l) 

+^ii((o)(e(p)i (g) (e(p)i +^i2(a))((e(p)i {^^)2 

+ (e(p)2^(e(p)i)+jf22(a))(e(p)20(e(p)2] (E), (22) 



where 



Xij=PiXij (/,7 = 1,2). 



(23) 



Notice that Xn = Xii^ so the susceptibility tensor x 
is symmetric, as would be required by the reciprocity 
considerations 
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FIG. 3: (Color online) (a) Real and (b) imaginary parts of the com- 
ponents of the effective dielectric permittivity tensor for the split- 
ring metamaterial with geometrical parameters given in the cap- 
tion of Fig. [21 In calculations, we neglect the force F'^^", so that 
k\2 = ^21 = and and ^22 are given by Eq. ( [T6l ). The copper ring 
is characterized by (Op = 2000 THz, 7=8 THz; permittivity of the 
dielectric is = 4. 



III. EFFECTIVE PARAMETERS OF A PCM 

A. Permittivity and permeability tensors 

To move on from a unit cell to the entire lattice of meta- 
atoms comprising a PCM, one needs to take into account that 
the field in each unit cell is modified by the presence of the 
neighboring meta- atoms. (It is known that an array of sym- 
metric meta- atoms with no intrinsic chirality can exhibit ex- 
trinsic chiral properties due to the way the atoms are arranged 
in a lattice.—) If the arrangement of meta-atoms is not too 
dense so that the individual atoms remain distinct, it can be 
assumed that the influence of the neighboring atoms is weak 
and can be simulated by regarding the meta-atoms as effective 
dipoles. The resulting field in each cell is equal to the sum of 
the electric field averaged over the whole planar metamaterial 
E and the field of the dipoles: (E) = E + 47rC(P), where C 
is the interaction matrix.— The interaction matrix depends on 
the symmetry of the lattice. For a planar arrangement of atoms 
(see Appendix A) its form is 



C = 0.36— (c;, 



) + 63; (g) - 2ez (8) e^) . 



(24) 



For the case D<^L, the interaction matrix has negligible com- 
ponents, and the average field in the metamaterial nearly co- 
incides with that in a single meta-atom: (E) ^ E. In our case 
D/L ^ 0.1, therefore, the influence of the surrounding meta- 
atoms can be non-negligible. 

The resulting effective dielectric permittivity tensor of the 
PCM can be derived from equations (E) = E + 4;rC(P), 
CeffE = E + 4;r(P), and (P) = ;t:(E). Its final form is 



£^ff{co) = l^47i{x ^-47tCy 



and can be rewritten as a matrix: 



^eff - 





^xy 





^xy 


£yy 














(25) 



(26) 



7 



The components £ij are complex, and it is seen that Cgff has 
the structure of a dichroic and anisotropic medium, in Hne 
with crystallographic expectations.^^ 

The typical frequency dependencies of e^x, %, and % for 
the example CSR structure in Fig. [2] are shown in Fig. [S] It 
can be seen that the structure features a Lorentz-like absorp- 
tion resonance in the range near 5.1 GHz, in agreement with 
experimental results for such CSRs.^^ This resonance corre- 
sponds to the minimum in the denominator Dq in Eq. ([T9b , so 
it is an intrinsic excitation in an individual meta-atom. 

It is important to realize two fundamental limitations of the 
presented homogenization approach. First, we neglect the 
magnetic dipole and electric quadrupole contributions (they 
have the same order of magnitude and should be accounted for 
simultaneously^^). This can be safely done for the light nor- 
mally incident onto a planar structure. Indeed, the quadrupole 
moment Q has the form Q{z) = (2j_ + Qzz^z ^ ^z) exp(i^o^). 
where Q^e^ = ^zQ± =0 and the z-dependence in Q{z) is 
caused by the external field. Hence the quadrupole contri- 
bution has the form V2(z) = ikoezQ{z) and is z-directed, re- 
sulting in a contribution to e^z- The magnetic dipole mo- 
ment is z-directed, too, since the electrons move in the (x, y)- 
plane, which results in effective magnetic permeability fl = 
diag(l, l,/!^^). Therefore, these higher-moment contributions 
will not affect £xx, £yy, and e^y, which are the only components 
that will play a part in determining the normal-incidence trans- 
mission and reflection spectra. So it is sufficient to consider 
just the electric dipole moment, assuming /i = 1 for the PCM. 

Secondly, and perhaps more seriously, the presented 
approach is commonly employed for bulk metamaterial 
homogenization,^^ and its applicability for metamaterial sur- 
faces leaves room for uncertainty regarding how, specifically, 
the homogenization in the z-direction should be performed. 
It is questionable whether a planar surface can be described 
as a finite-thickness slab of a bulk effective medium that 
would mimic the response of the metamaterial for all cases 
of the incident light, even when magnetic dipole and electric 
quadrupole contributions are taken into account. It is com- 
monly assumed that first-principle characterization methods 
based on multipole expansioni^ai^i^l should be used instead 
of homogenization. Still, we can obtain the effective material 
parameters valid in a specific case, in order to see whether a 
bulk material with planar chiral properties can be related to 
real PCMs.^^ Hence we can continue with the effective per- 
mittivity derived in Eq. (l26l) , keeping in mind that it is only 
valid for normal incidence of light. 



B. Polarization eigenstates 
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FIG. 4: Coefficient rf vs. frequency calculated for parameters indi- 
cated in Fig. [3] 



where n is the effective refractive index for the eigenwave. 
The electric fields in two eigenstates take the form 



E±=£x 1,- 



2£ 



. (28) 



Since the components of permittivity tensor are complex, 
these eigenwaves are elliptically polarized. Their direction of 
rotation ("right- or left-handedness") can be defined by the 
parameter t]± = {\£xy/Ex\)^iez{E± x E5_)3 An RH-polarized 
wave has ri± > 0, an LH-polarized one has r]± < 0, and, obvi- 
ously, T7± = coresponds to a linearly polarized wave whose 
sense of handedness cannot be determined. 

Consequently, we can identify whether the calculated effec- 
tive parameters of the CSR metamaterial correspond to those 
of a PCM by simply evaluating the product rf = t]^?]-. Con- 
ventionally, in isotropic or lossless birefringent media 7] = 
as the eigenwaves are linearly polarized. In 3D chiral or Fara- 
day media, i] < as the eigenwaves (either circularly or ellip- 
tically polarized, depending on the presence of anisotropy in 
addition to optical activity) are counter-rotating. On the con- 
trary, PCMs (and elliptically dichroic crystals, see Ref. 113) 
are characterized by co-rotating polarization eigenstates, so it 
is expected that T] > in these media. 

Fig.lHshows the coefficient T] calculated for the CSR meta- 
material in Fig. [21 with dielectric permittivity shown in Fig. [51 
Indeed, it can be seen that rf is positive for all frequencies in 
the vicinity of the intrinsic resonance around 5.1 GHz. Thus, 
the signature crystallographic property of a PCM (co-rotating 
elliptical polarization eigenstates) is indeed reproduced in the 
effective medium, confirming that it is 2D rather than 3D chi- 
rality that manifests in CSR metamaterials. 



Polarization of the eigenwaves of anisotropic medium with 
dielectric permittivity ([26l) can be found from the wave 
equation^ 



(^^(l-e,0e,)-eeff)E = O, 



(27) 



C. Transmission and reflection spectra 

As the final step in the theoretical model, we briefly out- 
line the calculation procedure for the optical spectra of a 
PCM. Following the set-up in Ref. 12, we consider the CSR 
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metamaterial of effective thickness D-\-h ^ D and calculate 
the reflection and transmission coefficients of a dichroic and 
anisotropic monolayer with dielectric permittivity given by 
Eq. (l25l) . Since we are interested in all possible polariza- 
tions of the incident wave, it is convenient to make use of 
the well known covariant operator generalization of the trans- 
fer matrix method (the covariant Fedorov's approach^). For 
the details on this method, the reader is referred to previous 
publications 

We define a unit vector q = pointing in the propagation 
direction and write the Maxwell equations for a monochro- 
matic normally incident wave in the form 



-i^o^E, 



q^-E = i/:oH, (29) 



where ko = CO /c is the vacuum wavenumber, and denotes 
the antisymmetric tensor dual to the vector q [(q^ )ii = (^ijkQj^ 
S'ijk is the Levi-Civita pseudotensor].^^ 

The fields are always tangential, and the field vectors are 
continuous across the layer interfaces. Eqs. (l29l) can be com- 
bined into the form 



dW(z) 
dz 



--ikoMW{z), 



where, for non-magnetic, non-gyrotropic materials. 



W = 



H 

qxE 



M = 







(30) 



(31) 



Here / = 1 — q(g)q = — q^ is the projection operator onto 
the plane normal to q, and 1 is the three-dimensional iden- 
tity tensor. The fundamental solution of Eq. (l3Ql) is a matrix 
exponential 

W(z) = P{z)W{0), P{z) = exp(i/:oMz), (32) 

where 4x1 dimensional constant vector W(0) is the initial 
field. The 4x4 matrix P{z) is called the evolution operator. 

Taking into account that in the medium surrounding the 
metamaterial layer (i.e., in air) the fields are related as q x E = 
±/H depending on the propagation direction,^ ^ the incident 
and reflected waves at the input (air/PCM) boundary are re- 
lated to the initial field W{0) as 



W(0) = 



Hi, 



ft 



refl- 



(33) 



From Eq. ([32b , the field at the output (PCM/air) interface is 
P{D)W{0). The evolution operator of the metamaterial P can 
be derived by setting £ = Cgff in Eq. dSTI) . 

On the other hand, the field at the output interface is the 
transmitted wave 



W(Z)) 



Htr- 



Hence the boundary conditions take the form 

/ 



Hi, 



Hrefl 



(34) 



(35) 



Multiplying Eq. ([35]) by the rectangular block matrix 
( / / ) and thus eliminating Hrefi, the expression for the 
transmitted magnetic field becomes 



Htr = 2 



{I I)p- 



Hi, 



(36) 



which, along with the evolution operator P, will depend on CO. 
Finally we define the transmission coefficient of the metama- 
terial slab as the ratio between the intensity of transmitted and 
incident waves: 



nco): 



|Htr(C0)P 
iHincP • 



(37) 



Eqs. ([36b and ([37b hold likewise for the electric fields.— 



IV. COMPARISON WITH NUMERICAL SIMULATIONS 

Varying the frequency and polarization of the incident wave 
Hinc , one can obtain the corresponding transmission spectrum 
as outlined in the previous section. We will be particularly in- 
terested in investigating the PCM transmittance for LH vs. RH 
circularly polarized incident wave, labelled Tl{co) and Tr{co), 
and corresponding to complex vectors Hinc = '^i^x iie^;), 
respectively. Here, Tl{co) and Tr{co) are the overall trans- 
mittances, without regard for polarization of transmitted light. 
The vast majority of materials (either naturally occurring or 
artificial) do not discriminate between LH and RH circular po- 
larization in transmittance, so that AT{co) = Tl{co) — Tr{co) is 
zero for all frequencies. A non-zero AT signifies the presence 
of circular dichroism and enantiomeric asymmetry. 

To test the applicability limits of the proposed model, we 
first compare analytically and numerically calculated Tl{co) 
and Tr{co) for our example CSR structure of Fig. [2] with vary- 
ing ring radius R and lattice period L. Numerical results 
are obtained using a commercially available finite integration 
solver (CST Microwave Studio) in the frequency domain, us- 
ing periodic boundary conditions in thcx — y directions. 

We see that a resonant dip in the transmission that results 
from the intrinsic resonance for the components of £ in the an- 
alytical model (see Fig. [3]) is reproduced in numerical calcula- 
tions and corresponds to the fundamental dipole excitation of 
the CSR. The resonant frequency /res changes when the ring 
radius is varied [Fig. Oa)]. The resonance has a Fano-like 
shape, which is also reproduced numerically. 

For frequencies below /res, we see a good agreement be- 
tween analytical and numerical results, which gradually wors- 
ens as R is increased in comparison with L, so that meta- 
atoms become closer to each other and the assumptions about 
a sparsely packed lattice that were needed in deriving Eq. ([25b 
become increasingly violated. This also causes a mismatch 
between analytically and numerically derived /res- On the 
other hand, for frequencies above /res the agreement is worse 
because the numerical spectra are affected by higher-order 
CSR resonances (which are explicitly not accounted for in our 
determination of a CSR's response), as well as the Bragg res- 
onances of the lattice, which are also neglected in our account 
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FIG. 5: (Color online) Comparison between the proposed model and the numerical simulation (CST Microwave Studio) for the transmission 
spectra for left-handed (Ti) and right-handed (Tr) circularly polarized incident waves. Calculations are done for asymmetric ring metamaterial 
with = 140°, a2 = 160°, /3i = 40°, and P2 = 20° for different ring radiuses and lattice constants: (a) for L = 15 mm and varying R, (b) for 
R = 6 mm and varying L. Other parameters are as given in Figs. [3 and [3] 
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FIG. 6: Same as Fig.[5]but for the transmission difference of AT = Ti — Tr. 



under the assumption that the effective medium is regarded as 
homogeneous. 

Fig. Ob) shows the dependence of the analytical vs. nu- 
merical spectra as L is varied for the constant GSR radius. 
It is seen that the agreement below /res is restored as L in- 
creases, confirming our reasoning. However, the agreement 
above /res becomes much worse because the Bragg resonances 
scale as /Bragg ~ 1/^, and are thus pushed into lower frequen- 
cies. Physically it means that the upper frequency limit where 
L <C A (and where the structure can be regarded as a meta- 
material) becomes smaller. Ultimately /Bragg moves past /res , 



which is where the response of an individual meta-atom be- 
comes irrelevant to the whole optical properties of the struc- 
ture. 

Therefore as far as the quantitative agreement of T{(o) is 
concerned, the model is found to be valid for the lattice pe- 
riod L not exceeding 20-25 mm and for the CSR radius R sig- 
nificantly lower than L. Most discrepancies occur in higher 
frequencies (/ > /res) and can be attributed to higher-order 
and/or lattice resonances which have been left out of consider- 
ation intentionally. The trade-off between split-ring coupling 
and grating diffraction is mentioned in the recent Ref . |41|. 



10 



Good agreement within these vaUdity Hmits of the model 
can also be confirmed in the spectra for the asymmetry 
AT{co) = Tl — Tr (Fig.[6l). Moreover, we see that in a vast ma- 
jority of cases, Tl ^ Tr everywhere except the vicinity of /res, 
which coincides with the range where polarization eigenstates 
are co-rotated elliptical (see Fig. [4]). So, many quantitative dis- 
crepancies in the transmission spectra have no effect over Ar 
and the model remains qualitatively valid for all the param- 
eter values shown in Fig. [6] with mismatch to the resonance 
frequency /res and the maximum value of Ar(/res) gradually 
increasing as the approximations behind the presented model 
become less accurate. The exception is the case of larger L 
where non-zero AT is also seen at odd-numbered higher-order 
CSR resonances [Fig. Ob)]. However, since such additional 
chiral response is spectrally well separated from the funda- 
mental resonance that interests us, it does not affect planar 
chiral properties of CSRs under present investigation. 

Hence, numerical simulations confirm that the proposed 
microscopic description of the CSR structures reproduces 
the PCM behavior, as reported in previous experiments 
Within the assumptions of the model that takes into account 
only the fundamental particle plasmon resonance of the CSR 
segments, the model provides a good agreement in a wide 
range of parameters. Having established this, we move on to 
investigate how the spectrum AT{co) behaves in various CSR 
designs. 

V. GEOMETRICAL TRANSFORMATIONS WITH CHIRAL 
SPLIT RINGS 

A non-zero Ar, signifying the presence of circular dichro- 
ism, carries special significance for 2D structures. One notices 
that spatial inversion of the whole system with respect to the 
plane normal to the metamaterial changes the handedness of 
the circular polarization (LH o RH), and replaces the struc- 
ture with its enantiomeric counterpart (see Fig.[T]). Hence, if 
T{(o) is the spectrum of any planar structure and T{co) is the 
spectrum of its enantiomeric counterpart, then 

TlA^) = (38) 

for any structure at any frequency. Therefore, for any planar 
meta-atom with no distinct 2D enantiomers (i.e., with an in- 
plane symmetry axis) Ar = 0. So, it is important to point out 
that Ar 7^ indicates the presence of planar chirality. A CSR 
of the design considered here (see Fig. [J]) becomes symmetric 
and therefore achiral if either the ring segments are of equal 
length (ai = a2), or the gaps between the segments are equal 
(j3i = ft), or else in a few degenerate cases when there is 
effectively just one segment (i.e., ai = 0, a2 = 0, j3i = 0, or 
i32 = 0). 

Furthermore, the reversal of the direction of incidence also 
transforms the structure into its enantiomeric counterpart but 
does not change the incident wave polarization. Hence, if the 
transmission spectra for the forward- vs. backward-incident 
wave are labeled T^'^, respectively, then Eq. (1381) results in 

T[i{0^) = fiiico) = T^iico), (39) 
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FIG. 7: (Color online) Difference of transmission for LH and RH 
circularly polarized incident waves for different ring geometries with 
L=\5 mm, R = 6 mm, oci = 140°, 0^2 = 160°, pi = 60° - fc: (a) 
comparison of AT for enantiomeric counterparts (pi ^ P2) in the 
range 10° < J82 < 50°; (b) evolution of maximum AT (Arpgak) be- 
tween symmetric and asymmetric CSRs for p2 < 30°. 



again, for any structure at any frequency i*^ Therefore, a non- 
zero Ar is a measure of the planar structure's directional 
asymmetry and its magnitude can be used as an indication of 
how strongly the planar chiral properties of a structure mani- 
fest themselves optically. 

Consider first the displacement of one of the ring segments 
along the circle by varying j3i and setting P2 = 360° — Pi — 
(ai + a2). The results are presented in Fig. [71 Not surpris- 
ingly, chiral properties are rather weak for Pi P2 and be- 
come larger as P2 increases towards the case of Fig. [5] where 
Pi = 40°, P2 = 20°. After a certain optimum value, however, 
Ar decreases again, vanishing completely in the symmetric 
case Pi= P2 = 30°. So the proposed theory confirms that both 
circular dichroism and directional asymmetry indeed vanish 
when mirror symmetry is present. 

The spectral shape AT{(o) is seen to have a wider shape 
for small P2, becoming the narrowest for the optimal case 
and then diminishing without significantly changing its shape. 
This is what one would expect as the inter- segment coupling 
(which is stronger for smaller P2 because the tips are in close 
proximity) pushes the particle plasmon resonances of a CSR 
apart from each other. This feature is specific to CSR de- 
sign: while the response of each arc-shaped segment in a CSR 
closely resembles that of a rod of equal lengthr^ the split 
ring is a geometry where the tips of the segments are in much 
closer proximity than for the rods placed at similar distance. 
Hence, the field enhancement near the tips causes the response 
of the whole CSR to depend strongly both on the individual 
segments and on inter- segment coupling. 

In addition, making P2 small while maintaining the CSR 
orientation should increase the contribution of extrinsic effects 
to chiral properties 1^ This is indeed seen in Fig. [3b), and this 
is likely the reason of a small but non- vanishing AT for these 
values. 

Displacing the ring segment past the symmetric case Pi = 
P2 = 30°, we notice in Fig.[7Ia) that Ar changes its sign and 
that an exchange of Pi ^ P2 results in the inversion of pla- 
nar chiral properties (AT{co) ^ —AT{co)). This entirely con- 
firms the result expected from Eq. (1381) , taking into account 
that structures obtained by an exchange of Pi and P2 are 2D 
enantiomeric counterparts. 

We also examine CSRs with variable length of the longer 
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FIG. 8: (Color online) Difference of transmission for LH and RH 
circularly polarized incident waves: (a) for ai = 140°, p2 = 20°, 
varying pi and 0^2 = 200° - Pi ; (b) for 0^2 = 160°, = 20°, varying 
pi and ai = 180° 



ring segment by varying Pi and a2 while keeping the other 
two angles constant (note that all the time ai + a2 + j3i + j32 = 
360°). The results are given in Fig.Oa). As expected, Ar van- 
ishes for the two symmetric cases j3i = ft = 20° and Pi = 60° 
(ai = a2 = 140°). The sign of Ar changes when these two 
symmetric cases are traversed. The increase of the peak fre- 
quency /res corresponds to an increase in the stiffness coeffi- 
cients ^11^2 as one of the segments becomes shorter, affected 
by a change in the coupling between the segments' resonances 
as they differ in length more strongly. 

Similarly varying the length of the shorter ring segment, 
i.e., changing the angles Pi and ai, the same behavior is ob- 
served, as can be seen in Fig.[8tb). The spectral shape changes 
in the same manner as in the previous case, the peak in Ar be- 
coming broader for larger Pi . However, this case is specific 
because 0^2 + & = 180°, so there is only one symmetric shape 
(Pi = & = 20°, ai = a2 = 160°). Because of this "degener- 
acy", Ar does not change sign when traversing the symmetric 
case. The same property can be responsible for higher peak 
values of Ar(/res). 

Finally, to determine how the chiral properties of PCMs 
scale with the size of the meta-atoms, we investigate the de- 
pendence of Ar on the ring radius R in Fig. |9l The resonant 
frequency /res depends on the radius primarily due to 1/7? ap- 
pearing in the stiffness coefficients kij [see Eq. ([T6l) 1 . The es- 
timated dependence /res ~ l/Ris confirmed in Fig. Ob). The 
resonance peak in AT{co) broadens as R increases. 
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FIG. 9: (Color online) (a) Difference of transmissions of LCP and 
RCP incident waves for different values of ring radius R. (b) Depen- 
dence of the peak frequency /res (7?). 



VI. CONCLUDING REMARKS 

In conclusion, we have proposed a microscopic theoreti- 
cal description of planar chiral metamaterials based on the 
electronic Lorentz theory. Using a chiral split-ring (CSR) 
geometry as an example and considering the dynamics of 
individual free electrons, we arrive at expressions for elec- 
tromagnetic response of a single split-ring meta-atom. Its 
polarizability is derived analytically without phenomenolog- 
ical parameters. The effective dielectric permittivity tensor 
£eff is then obtained from the single-atom response along the 
lines of standard homogenization techniques, -^2^ to the ex- 
tent that these techniques can be used at all to describe planar 
structures.—"— Even in this simplified consideration, the crys- 
tallographic structure of this tensor is shown to coincide with 
that of elliptically dichroic media, as would be expected from 
earlier theoretical studies.— 

The transmission spectra of a GSR-based PGM are then cal- 
culated using the standard wave operator based extension of 
transfer-matrix methods. ^^'^^ In agreement with previous the- 
oretical and experimental results, ^^'^^ the spectra show a dif- 
ference with respect to whether the incident wave has left- or 
right-handed circular polarization [AT{co) = Tl{co) — Tr{co)]. 
This difference, which translates to directional asymmetry, 
is shown to strongly depend on the geometrical parameters 
of the GSR (Figs. |7]-[8]). Whenever the split ring is symmet- 
ric, i.e., when there are no distinguishable 2D enantiomers, 
optical manifestation of planar chirality is seen to vanish 
[AT{co) = 0], and Ar changes its sign when the structure is 
replaced with its enantiomeric counterpart [Fig.Oa)]. 

Note that we have deliberately chosen the overall GSR ori- 
entation in the lattice so as to focus on purely intrinsic chirality 
and to suppress extrinsic effects where possible. A detailed in- 
vestigation of how intrinsic and extrinsic chiral properties in- 
teract in GSR-based PGMs warrants a separate investigation. 

While the specific split-ring geometry is chosen for its rel- 
ative simplicity in analytical derivation, it should be under- 
stood that the proposed approach can be extended to any pla- 
nar meta-atom consisting of thin wire-like metallic elements 
where transverse motion of electrons is restricted. The inte- 
grals in Eqs. ^ and (fTTI) are likely to be more complicated 
and may have to be taken numerically. Moreover, it may be 
particularly challenging to determine the correct charge den- 
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sity dependencies [see Eq. (|2])] and to identify the loops that 
contribute to the effective mass; one may need to use equiv- 
alent LC-circuit or even resort to using data from direct nu- 
merical simulations if the geometry is particularly intricate. 
Nevertheless, once the equations for the effective mass and 
stiffness coefficients for a particular unit cell geometry are es- 
tablished, parametric transformations of this geometry lend 
themselves to very easy semi-analytical treatment within the 
proposed framework. 

Furthermore, oblique wave incidence and non-planar meta- 
atoms composed of similar thin elements are also tractable if 
the magnetic dipole and electric quadrupole contributions are 
accounted for, giving rise to corrections in Cgff and introduc- 
ing effective magnetic permeability tensor /i, as well as gyra- 
tion pseudotensors responsible for magnetoelectric coupling 
or spatial dispersion i^Si^Lffi However, it still remains an open 
question whether such oblique-incidence treatment would be 
universal taking into account the inherent limitations of apply- 
ing volume homogenization to surface structures.—"— Should 
such a generalization prove feasible, it is very interesting to 
extend the proposed approach from a single-layer PCM to 
PCM-based multilayers (in the cases when such PCMs can 
lend themselves to 3D homogenization^^) and investigate its 
applicability as a planar metamaterial turns into a bulk one. 
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field in an ^th unit cell as a sum of the averaged electric field 
in the metamaterial E and the fields of the dipoles p at the 
center of the each cell: 



(E)=E- 



I 



3ri 



-r? 



-P, 



(Al) 



where is the radius-vector of /th dipole. Following Ref. l2lL 
we define the interaction matrix C by means of equation 
(E) = E + 47rCPo- The polarization of the medium is con- 
nected with the dipole moment of the single cell via Pq = A/op, 
where A/q = 1 / {L^D) is the number of inclusions per unit vol- 
ume (we suppose that the cell is square). Therefore, the inter- 
action matrix equals 



(A2) 



Placing the origin at the center of the cell under consider- 
ation (^th cell) we present the radius-vectors of the others as 

= ^kl = L{kex + /e^), where integer numbers k, I vary from 
— oo to oo. The nth cell is characterized by the numbers ^ = 
and / = and should be excluded from the summing. 

The result of the summation is used as Eq. (l24l) in Sec. Hill 
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C= — (/-2e^(g)e^), 



(A3) 



where / = + e^; is the projection operator onto the 
plane normal to the z axis. 



Appendix A: Derivation of the interaction matrix 



2n 



: 0.3592 (A4) 



To calculate the interaction matrix for meta-atoms arranged 
in a square planar array used in Sec. IIII Al one can write the 



and ^(x) = ^ ^ is the Riemann Zeta function. 
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